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i Abstract. The characteristic formulation of the relativistic hydrodynamic equations 

^ \ (Donat et al 1998 J. Comput. Phys. 146 58), which has been implemented in many 

relativistic hydro-codes that make use of Godunov-type methods, has to be slightly 

. modified in the case of evolving barotropic flows. For a barotropic equation of state, 

| a removable singularity appears in one of the eigenvectors. The singularity can be 

£Nl ' avoided by means of a simple renormalization which makes the system of eigenvectors 

£^ , well defined and complete. An alternative strategy for the particular case of barotropic 

' flows is discussed. 

o 

(N 1 PACS numbers: 04.25.D-, 47.11.-j, 47.75. +f, 95.30.Sf 



On Numerical relativistic hydrodynamics and barotropic equations of state 



2 



Barotropic fluids obeying an equation of state (EoS) of the type p = p(e) (barotropic 
EoS), where e = p(l + e) is the total mass-energy density (p and e being, respectively, 
the rest-mass density and the internal specific energy), are a particular class of fluids 
which become familiar, for example, to cosmologists interested in the analysis of linear 
perturbations of radiation-dominated universes (see, e.g., [2]), or to astrophysicists 
interested in the physics of dense matter where the asymptotic EoS for high density 
and cold matter is barotropic (see, e.g., [3]). The mathematical properties of relativistic 
Euler equations for barotropic EoS have been studied in [1]. A particular case of a 
barotropic EoS is the ultrarelativistic EoS, where pressure is directly proportional to 
the energy density. Recently, several authors have paid particular attention to this 
ultrarelativistic EoS in order to find the exact solution of the relativistic hydrodynamical 
Riemann problem with nonzero tangential velocities [5], the corrugation instabilities 
occurring for solutions of that Riemann problem [6], or just in order to compare 
efficiencies of relativistic hydro-codes specifically designed for such particular EoS [TJ. 

Twenty years ago, authors in [8] paved the way for modern numerical relativistic 
hydrodynamics (RHD) extending to this field the strategy of using high-resolution shock- 
capturing (HRSC) methods. By exploiting the hyperbolic and conservative characters 
of the RHD equations, these authors (and collaborators) derived the full spectral 
decomposition of the Jacobian matrices associated with the fluxes in a series of papers: 
[8], in ID (one-dimensional) general relativity hydrodynamics (GRHD); [9] and [T], in 
multidimensional special relativity hydrodynamics (SRHD); [10] (where the eigenvectors 
were given for a diagonal metric) and [H] (where authors derived the full spectral 
decomposition: right and left eigenvectors) in multidimensional GRHD. This numerical 
strategy has proved to be very fruitful in the field of computational RHD (see, e.g., the 
reviews [12] and [13]). 

All the above papers (and works related with them) as well as all the current 
numerical codes based on these references make use of a general EoS of the form 
p = p(p,e). This general EoS contains some particular cases (see, e.g., [TJ], [15], [16]): 
(i) dust (p = 0). For this EoS, motion equations can be reduced to the inviscid Burger's 
equation (in ID), (ii) A gamma-law EoS (p = Kp 1 ). This case makes one of the motion 
equations (the energy equation) redundant. 

This note is addressed to those researchers, working in relativistic astrophysics 
or numerical relativity, that might be interested in studying the evolution of barotropic 
fluids by using relativistic hydro-codes originally designed to use Godunov-type methods 
and a general EoS, of the type p = p(p, e). 

The evolution of a relativistic fluid is described by a system of equations which are 
the expression of local conservation laws: the local conservation of baryon number and 
energy-momentum, respectively, 



V„J" = 0, 
V^T"" = 0, 



(1) 
(2) 
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the current J M and the energy-momentum tensor T^ u being 

r = pu», (3) 
= phu^Uv + pg^. (4) 

(Greek and Latin indices run, respectively, from to 3 and from 1 to 3. In the units we 
use, the speed of light is equal to one). 

In the above equations, h is the specific enthalpy, defined by h = 1 + e + p/p, is 
the 4-velocity of the fluid, and g^ v defines the metric of a general spacetime M. where 
the fluid evolves, V M stands for the covariant derivative. 

In the following we will restrict our analysis to the Minkowski spacetime. The above 
system of equations, ([I]) and (J2J), can be written as: 

dF°(w) aF'(w) _ dF°(w) pt «9F°(w) _ 

where 

9F*(w) 

are three 5x5 Jacobian matrices and w = (p, e) is the 5- vector of primitive variables. 

v? 

The components of the 3- velocity, v 3 \ are defined according to v 3 = and the Lorentz 

u° 

factor, defined by W = it , satisfies the familiar relation = (1 — v 2 ) -1 / 2 , where 
v 2 = r]ijV % v 3 \ with if 3 being the spatial metric of Minkowski spacetime. 
The components of the 5- vector of conserved variables, 

F°(w) = (D, S 3 , r) = (pW, phW 2 v j , phW 2 -p- pW) , (7) 

are, respectively, the rest-mass density, the relativistic momentum density and the total 
energy density without the rest-mass density. 
The three 5-vectors of fluxes are: 

F(w) = (pWv\ phW^tfv* +pv ij , phW 2 v l - pWv 1 ) . (8) 

Assuming an EoS of the form p = p(p, e), the local sound velocity c s is given by 

^ = X+(^)« > X=| • ( 9 ) 

Authors in [1] showed the characteristic structure of the multidimensional RHD system 

of equations. Some components of one righteigenvector, and also the determinant of 

the matrix of the righteigenvectors (Eqs. (17) and (21), respectively, in pQ) contain 

t\j 1 dp 

the expression K := with k := — — . It is clear that K is not defined when 

k — c s p oe 

k — c 2 = 0, and therefore, in this case, the system of righteigenvectors given in [T] is not 
well behaved. Regarding this point, the following proposition can be proved. 
Proposition The thermodynamic function k — c 2 is identically zero iff p — p (e), e := 
p(l + e), with p being an arbitrary function. 
Proof. 
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k — c s = 0, taking into account the definition of c 2 , can be shown easily to be 

equivalent to (1 + e) k — x — 0. Introducing the auxiliary variables t := log(l + e) 

r, Op dp 

and x := logo, we conclude that k — c 2 = — — — — = 0. The general 

at Ox 

solution of this advective-like partial differential equation (PDE) is p = p D (x + t), with 
x + t = log(p(l + e)) = loge and being p Q an arbitrary function. 

This conclusion can be also obtained in terms of functionals. k — c 2 = 0, from the 

3c Op dc dp 

definitions of c 2 and e, is equivalent to — — -rr-^r = 0, which describes a functional 

Oe Op Op Oe 

dependence of p and e, p = p (e). 

Hence, we have showed that k — c 2 = ^ p = p (e) . □ 

From the above proposition, we conclude that the characteristic eigenvectors of the 
SRHD equations (J5]) derived in [1] have to be reviewed just in the case of barotropic 
EoS. 

For the sake of conciseness, let us concentrate on the ID version of the analysis 
presented in [T]. The expressions for the righteigenvectors are: 
/ K _ K\ 

r = 7777,M- 7777 , (10) 



hW 1 ' hW t 

r ± = (1, hW(v±c s ), hW(l±vc s ) - 1) . (11) 

Multiplying ro by the normalization factor k — c 2 s , the determinant of the matrix of 
renormalized righteigenvectors, Aid, satisfies 

A 1D = %^. (12) 

The singularity has been removed, both in the eigenvectors and the determinant, and 
the system of eigenvectors is complete for a general EoS, including a barotropic one. 

Analogously, in the multidimensional case, it can be showed that, given an 
arbitrary spacelike 3-vector Q, the general eigenvalue problem for the matrix QB^ (see 
equation (jSJ)) has the following renormalized righteigenvectors: 

^=(^, -A), ( 13 ) 

102,03 = {W Vj ( j , h(C + 2W\( j v i ), (2hW - l)WvjC j ) , (14) 
r± = (l, hW (^ A -^ l) C + v^j , hWA ± - lj , (15) 

where A± = — — ; tt 7 , and denotes the two unitary spacelike 3- vectors which 

together with form an orthonormal basis. The particular values of the above 

righteigenvectors in the direction Q = (1,0,0) coincide, as it should be, with the 
ones derived in pp. The singularity has been removed and all the components of the 
eigenvectors are well defined. The determinant of the righteigenvectors, A, does not 
vanish, 

-clWh\\ + -\^{{v^f-^f 
(A. r+J c ; V)(A r^-QQ)^ ' 
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and the system of eigenvectors is complete (even for a barotropic EoS). 

Therefore, the main conclusion is that, by an appropriate renormalization of one of 
the eigenvectors derived in [I], it is possible to use the corresponding numerical hydro- 
code for a general EoS, including the case of barotropic EoS. Similar renormalization 
is required in RMHD when the characteristic eigenvalues of the system of equations 
become degenerate [T7] . 

A second strategy is possible: to rewrite down the SRHD equations for a barotropic 
EoS. It is the one proposed in [5j El Ej. Yes, indeed, quoting Anile in [16] (page 11): 
'Since the rest-mass density does not intervene in the state equation, these fluids can be 
described solely by the conservation of energy and momentum equations. The mass 
conservation equation decouples from the others and can be solved after the fluid 
motion has been determined'. The above Anile's quotation defines the (alternative) 
computational strategy to follow. 

The conservation of momentum and energy equations governing the evolution of a 
barotropic fluid, i.e. the one satisfying a barotropic causal EoS, are 

d t [(e + p)W 2 v 1 } + ^[(e + p)W 2 v i v j + pr] ij ] = 0, (17) 
d t [(e + p)W 2 - p] + dj[(e + p)W 2 v j ] = 0. (18) 

The new vector of conserved variables is U = ((e + p)W 2 v l , (e +p)W 2 —p). The system 
of equations ( I17j) and ( I18j) can be cast into the form of a quasi-linear system of PDEs: 



d t U + A j djU = 0, A j 





Ai 


Si 






(19) 



where A{ = v J ' 5\ + [(1 — v 2 c 2 ) b~\v % + 2c 2 v j v k v l — 2c 2 v k rf j ] (1 — v 2 c 2 ) 1 and A j 2 = 
[-(l + c 2 s )v j v i + c 2 (l+v 2 )ri l ^}(l-v 2 c 2 s y 1 . 

Let (j be an arbitrary spacelike vector. The eigenvalues of the matrix QA? are 

Ao = Qv j , (20) 
x (1 - ± Cs(l - v 2 ) 1/2 ((c 2 - lXQV) 2 + (1 - rVR A T 2 , 91 , 

1 — v z Cg 

where Ao has multiplicity 2. The corresponding eigenvalues associated with previous set 
of eigenvalues are 

I > (22) 



r i (C^'-A ± )C + (A ± C^-0C> 4 i m] 



where Q denotes, as previously, the two unitary spacelike 3- vectors which together with 
form an orthonormal basis. The system of eigenvalues form a complete basis. 



Therefore, the system of equations ffl9l) is hyperbolic. 

In brief, according to this second strategy and for a barotropic EoS, the numerical 
procedure, regarding the use of Godunov-type methods, would consist in solving, first, 
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the system ( 1T9"j) . The rest-mass density can be obtained, in a second step, by solving the 
continuity equation in the field velocity provided by the first step. Quoting Anile in [16] 
(page 11): '...the only role of the continuity equation is to provide a volume tracer.' 

Summary 

The characteristic formulation of the RHD equations carried out in [1] needs a 
minor modification for a barotropic EoS in order to remove the singularity in one of the 
eigenvectors. 

Any user of a relativistic hydro-code, based on Godunov-type techniques and 
designed for a general EoS p = p(p,e), that might be interested in using a barotropic 
EoS, can follow one of the following two strategies. (A) Just to renormalize one of 
the eigenvectors, such as we have shown above. (B) To rewrite explicitly the SRHD 
equations for a barotropic EoS. 

In terms of efficiency, it is obvious that strategy (B) is the most recommendable, 
since there is one equation less in the system. However, we recommend strategy (A) due 
to the fact that it can be applied in studies of flows obeying a general EoS, including 
the barotropic ones. Moreover, strategy (A) allows one to handle with those flows that, 
during their evolution, might change their thermodynamical properties in such a way 
that the EoS admits a barotropic branch as, e.g., the asymptotic behavior of the EoS 
for dense matter at high densities. In such a case, strategy (A) would be the only one 
suitable for a global description of the evolution of these flows. 
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